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ABSTRACT 

Detailed spectral analysis of the Galactic X-ray background emission, or the Galactic Ridge X-ray 
Emission (GRXE), is presented. To study the origin of the emission, broad-band and high-quality 
GRXE spectra were produced from 18 pointing observations with Suzaku in the Galactic bulge region, 
with the total exposure of 1 Ms. The spectra were successfully fitted by a sum of two major spectral 
components; a spectral model of magnetic accreting white dwarfs with a mass 

of 0.6618$ Mq, and a 

softer optically-thin thermal emission with a plasma temperature of 1.2 — 1.5 keV which is attributable 
to coronal X-ray sources. Whe n combined with previou s studies which employed high spatial resolution 
of the Chandra satellite (e.g. Revnivtsev et al.]|2009[), the present spectroscopic result gives another 



strong support to a scenario that the GRXE is essentially an assembly of numerous discrete faint X-ray 
stars. The detected GRXE flux in the hard X-ray band was used to estimate the number density of 
the unresolved hard X-ray sources. When integrated over a luminosity range of ~ 10 30 — 10 34 erg s" 1 , 
the result is consistent with a value which was reported previously by directly resolving faint point 
sources. 

Subject headings: X-rays: diffuse background - Galaxy: bulge - novae, cataclysmic variables 



1. INTRODUCTION 



Tanuma et al. 1999 



Apparently extended X-ray emission has been observed 
along the Gal actic plane since early years of X-ray astro- 



physics (e.g. Cooke et al. 1969 Worrall et al. 1982) 



and its origin has been one ot the long-standing mys- 
teries in the research field. After its spatial structure, 
the emission has been called the Galactic Ridge X-ray 
Emission (GRXE). Its total luminosity is estimated to 
be ~ 1 — 2 x 10 38 erg s _1 in the conventional X-ray en- 
ergy range of 2 — 10 keV (e.g. Koyama et al.|1986 Valinia 



& Marshall 1998 ). The detection of intense emission lines 



from highly ionized Fe in the GRXE spectrum indicated 



Koyama et al. 


198G 


Yamauchi 


; was still unc 


ear whether the 



plasma is gravitationally bound to discrete sources, i.e. 
X-ray stars, or is of truly diffuse nature filling the inter- 
stellar space. 

Many scenarios which employ different origins have 
been proposed to explain the GRXE. They can be di- 
vided into two major groups depending on the assumed 
nature of the source, namely "Point Source" and "Dif- 
fuse" scenarios. The former assumes that a collection 
of faint but numerous discrete X-ray point sources in 



the G alaxy composes the GRXE (e.g. Revnivtsev et al. 



20061, similar to the case of the Cosmic .X.-ray Back 
ground (e.g. Giacconi et al. |1979| Ueda et al. 1999 
Mushotzky et al. 2000). In the latter scenario, m con- 
trast, literally-diffuse X-ray emitting materials are con- 
sidered to fill the interstellar space over the observed 
scales of the GRXE, i.e. is several tens of degrees (longi- 
tudinal) by a few degrees (latitudinal) around the Galac- 
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tic center (e.g. Kaneda et al. 1997 
and references therein ) . 

Unlike the "Diffuse" scenario which confronted several 
difficulties of energy supply and/or plasma confinement, 
almost only one uncertainty associated to the "Point 
Source" scenario was whether there are enough number 
of faint Galactic X-ray sources that can collectively ex- 
plain the GRXE surface brightness and spectrum, while 
individually satisfy the luminosity limit (< 10 33 erg s _1 ) 
not to be spatially resolved by X-ray telescopes. In previ- 
ous studies which used instruments with moderate angu- 
lar resolutions (e.g. > 1'), it was not possible to obtain a 
definitive conclusion on this problem from direct source 
counting. 

Chandra provided a great leap in the GRXE research 
thanks to its ever b est angular resolution in the X-ray 



wavelength (~ 0.5"). Revnivtsev et al. (2009) conducted 
a very deep Chandra observation on the Galactic bulge 
region (l,b) — (0.113,-1.424), and showed that more 
than 80% of the 0.5 - 7 keV GRXE surface brightness 
can be resolved into faint X-ray stars with a limiting 
sensitivity of ~ 10~ 16 erg cm~ 2 s _1 . Based on X-ray 
spectral hardness of t he faint sources , as we ll as their in- 
frared identifications, [Ebisawa et al. (2005) argued that 
most sources which show soft and hard spectral indices 
can be regarded as X-ray emitting stellar coronae and ac- 
creting white dwarf binaries like polars and intermediate 
polars (IPs), respectively. This has been also supported 
by a study of low-lu minosity X-ray sourc e population 
near the solar system ( |Sazonov~ et al. 2006 1_ 



The observation by |Revnivtsev et al. 1 2009 ) has in- 
vested ~ 1 Ms of Chandra exposure in a selected field. 
This is considered to be an ultimate limit achievable 
with the current state of the art. However, questions 
still remain: "how do individual types of X-ray sources 
contribute to the GRXE?" and "is there any additional 
spectral component in its spectrum that cannot be at- 
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tributed to these types of sources?" . 

To address these issues, we take an alternative ap- 
proach; broad-band and high-energy-resolution spec- 
troscopy of the GRXE. For understanding competitive 
contributions from multiple types of X-ray sources, it is 
essential to analyze the GRXE spectra on a broad-band 
basis because some of the putative constituents, namely 
magnetic accreting WDs, are known to emit strong sig- 
nals in the hard X-ray energy band above 10 keV. High 
energy resolution is also necessary to understand the 
characteristic Fe emission lines which are emitted pre- 
sumably by several types of sources. 

Some authors already examined broad-band spectral 
similarities between the GRXE and several types of 



X-ray sources (Revnivtsev et al. 2006 Krivonos et al 



2007 1 although spectral quality were not necessarily high 



Compared to these reports, our study improves data 
quality in terms of counting statist ics and energy resolu - 
tion by using the Suzaku satellite ( Mitsuda et al.| [2007). 
In the following sections, we try to spectroscopicaliy de- 
compose the GRXE along the following approach: 

(i) Construct a broad-band GRXE spectrum using the 
data of Suzaku Galactic bulge observations. 

(ii) Constrain contributions from magnetic CVs by an- 
alyzing the hard X-ray spectral shape of the GRXE 



et al. 


(2010 


) and 
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(2011 



(iii) Extrapolate the magnetic CV contribution down to 
a softer energy band below 10 keV, and quantify 
contribution from other types of X-ray sources to 
the soft-band GRXE. 

2. OBSERVATIONS AND PREPARATION OF DATA SET 

2.1. The Suzaku X-ray Observatory 

Suzaku, which is the fifth Japanese X-ray astrophysical 
observatory launche d in 2005, carries the X-ray Imaging 



Spectrometer (XIS; [Koyama et al. 2007) and the Hard 

XUT 



X-ray Detector (HXD; jl'akahashi et al.||2007[ ) 



The XIS consists of four X-ray imaging charge coupled 
device (CCD) camer as as a focal plane detect or of the X- 
ray Telescope (XRT; |Serlemitsos et al.|2007| ). The nom- 
inal energy covered by the XIS is 0.2 — 12 keV (Koyama 
et al.|2007 ), and the CCD provides a rectangular field of 
view (FOV) of 18' x 18'. 

The HXD is a non-imaging collimated X-ray detector 
which consists of stacked two main detection parts; sili- 
con p-intrinsic-n diodes (hereafter PIN) and gadolinium 
silicate (GSO) scintillation crystals. They cover energy 
ranges of 10 - 70 keV (PIN) and 40 - 600 keV (GSO). 
In the PIN energy band, a passive collimator provides 
an FOV of 34' x 34' at full-width at half maximum 
or a bottom -to-bottom aperture of 68' x 68' (Kokubun 
et al.||2007|). This tightly collimated FOV provides one 



of strong points of the HXD in studies of apparently ex- 
tended emission, because contamination from bright dis- 
crete X-ray sources is considerably lower than in previous 
hard X-ray instruments. 

Thanks to the low-earth orbit (altitude of ~ 570 km) of 
Suzaku, the XIS has lower non X-ray backgrounds (NXB) 
compared to XMM-Newton and Chandra, and is suitable 
for observing extended emission with moderately low sur- 
face brightness such as the GRXE. Simultaneous and 



broad-band spectral coverage below and above 10 keV, 
achieved with the two detectors, allows us to accurately 
measure spectral parameters of the GRXE, including in 
particular the Fe-K lines and hard X-ray continua, re- 
spectively. 



2.2. Mapping of the Galactic Bulge Region 

Since the launch, Suzaku has been actively perform- 
ing observations of the Galactic center and the Galactic 
bulge regions (\l\ < 5°,|6| < 3°) to study the extreme en- 
vironment near the supermassive black hole, hot plasmas 
around supernovae, molecular clouds reflecting X-rays, 
and a nu mber of X-ray point-like sources in these re - 
gions (ej;. Nobukawa et al. 2008 Uchiyama et al. 2011b ). 
Figurelljshow a mosaic image produced by collecting all 
these Suzaku data with a total exposure of 4.4 Ms. The 
data are taken from 92 mapping observations performed 
during 2005 — 2010 of which individual log is summarized 
for reference in Table |6] 



2.3. Selecting Suitable Observations 

To avoid signals from bright point sources which could 
contaminate GRXE signals, we carefully selected obser- 
vations that do not contain any bright poin t source. 



The INTE GRAL General Reference Catalog (Ebisawa 
et al. 120031) and the forth INTEGRAL/IBIS catalog ( [Hire? 



et al. 



2010p were used as references to positions and typ 



ical intensities of known X-ray sources. We discarded 
an observation when it is contaminated by X-ray point 
sources of which flux is above 0.1 mCrab (equivalently 
~2x 10 -12 erg cm~ 2 s _1 in the 2 — 8 keV) inside the 
HXD FOV. This criterion was set because a GRXE flux 
observed with the HXD is typically 1-2 mCrab in this 
sky region, and the limiting HXD/PIN flux of 0.1 mCrab 
corresponds to 5 — 10% of the detected GRXE flux. Ta- 
ble [I] lists observations which passed the criterion, and 
usedin the present spectral analysis. In Figure [l] (b), 
the FOVs of the utilized observations are overlaid on the 
same X-ray image as panel (a). To study possible spatial 
variations of spectral properties of the GRXE, we divided 
the accepted observations into two groups, and hereafter 
refer to them as Region 1 (on-plane, I < regions) and 
Region 2 (off-plane, b < regions). Compared to typical 
Suzaku observations, the total effective exposures of the 
selected regions are extremely deep, 592.1 ks (XIS) and 
524.0 ks (HXD) for Region 1, and 418.8 ks (XIS) and 
340.2 ks (HXD) for Region 2. 

2.4. Contamination from Detected Point Sources 

Although we selected observations with very low point- 
source contamination, we still notice four faint X-ray 
point sources, as indicated with green crosses in Figure 
[T] in the FOVs of Regionl. Source a is a newly- found 
soft X-ray source, and the remaining thr ee are previ- 
ously known sources; b=AX J17 42.6-3022 dSakano et~aL 
2002[) c=Suzaku J1740.5-3014 flUchiyama et al.||2011a| r 



aHdd=IGR J17391-3021 (c.g 
extracting the GRXE spectra 
sources are masked. On the other hand 



ozzo et al |2UlUp . When 
rom the XIS data, these 



we cannot ex- 
clude them in HXD/PIN because it does not have imag- 
ing capability. As a result, the HXD/PIN spectra contain 
signals from these point sources, in addition to those from 
the GRXE. Therefore, before actually extracting GRXE 
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TABLE 1 

Summary of the Galactic bulge observations. 



Obs. ID a Coordinate 11 Start time Exposure Count rate d 

I b UT XIS HXD XIS HXD 



Region 1 



1 


501053010 


-1.83 


-0.00 


2006-10-10 21:18:59 


21.9 


19.9 


0.25 


0.07 


2 


503014010 


-2.10 


-0.05 


2008-09-18 04:46:49 


55.4 


51.2 


0.18 


0.05 


3 


503015010 


-2.35 


-0.05 


2008-09-19 07:33:05 


56.8 


52.8 


0.19 


0.04 


4 


503016010 


-2.60 


-0.05 


2008-09-22 06:47:49 


52.2 


49.3 


0.18 


0.03 


5 


503017010 


-2.85 


-0.05 


2008-09-23 08:08:10 


51.3 


48.6 


0.17 


0.04 


6 


503021010 


-1.62 


0.20 


2008-10-04 03:44:03 


53.8 


49.6 


0.24 


0.07 


7 


503076010 


-1.50 


0.15 


2009-02-24 17:04:51 


52.9 


43.8 


0.27 


0.07 


8 


503077010 


-1.70 


0.14 


2009-02-26 01:01:00 


51.3 


43.7 


0.24 


0.07 


9 


504001010 


-1.47 


-0.26 


2010-02-26 09:15:00 


51.2 


42.2 


0.20 


0.05 


10 


504002010 


-1.53 


-0.58 


2010-02-27 16:14:41 


53.1 


46.6 


0.17 


0.04 


11 


504003010 


-1.45 


-0.87 


2010-02-25 04:33:17 


50.9 


41.3 


0.19 


0.02 


12 


504090010 


-1.49 


-1.18 


2009-10-13 04:17:20 


41.3 


35.0 


0.20 


0.03 










Total Exposure 


592.1 


524.0 






Region 2 


1 


502004010 


0.17 


-1.00 


2007-10-10 15:21:17 


19.9 


18.8 


0.45 


0.05 


2 


502059010 


-0.00 


-2.00 


2007-09-29 01:40:51 


136.8 


110.5 


0.35 


0.02 


3 


503081010 


0.03 


-1.66 


2009-03-09 15:41:50 


59.2 


57.6 


0.49 


0.01 


1 


504050010 


0.10 


-1.42 


2010-03-06 03:55:37 


100.4 


80.5 


0.60 


0.02 


5 


504088010 


-0.00 


-0.83 


2009-10-14 11:30:56 


47.2 


32.6 


0.43 


0.05 


6 


504089010 


-0.05 


-1.20 


2009-10-09 04:05:59 


55.3 


40.2 


0.54 


0.02 










Total Exposure 


418.8 


340.2 







a Observation ID. 

b Aim point in the Galactic coordinate (degree). 
c Net exposure in units of 10 s. 

d Count rates in units of counts s — 1 calculated over 1 — 9 and 17— 50 keV per one of the three 
XIS sensors and the HXD, respectively. NXB is subtracted. The XIS count rate includes 
CXB which accounts about 0.04 counts s — 1 whereas that for the HXD (~ 0.01 counts s — 
is subtracted in the HXD rates. 



spectra, we study their effects on our HXD/PIN GRXE 
data. 

We produced spectra of the four point sources from 
the XIS data of Obs.IDs 504001010, 503021010, and 
402066010, as shown in Figure [2] To estimate their 
fluxes in the HXD/PIN energy band, we fitted the spec- 
tra with phenomenological models. Because there is no 
characteristic feature in the spectra of Sources a,b, and 
d, we applied to them power-law models modified by 
the interstellar absorpti on (i.e. wabs xpowerlaw in the 
XSPEC fitting package; Arnaud 1996). Source c shows 
clear emission-line features at 6 — 7 keV, and is identi- 



fied as an intermediate polar by Uchiyama et al. (2011a I 
based on its periodic time variation. Therefore, this 
source should be fitted with a multi-temperature ther- 
mal plasma model. Because the counting statistics are 
not so high, we restricted ourselves to a model consist- 
ing of two collisional ionization equilibrium (CIE) plasma 
components plus a gaussian subject to photo absorptions 
[i.e. wabsx(APEC+Gaussian)+wabsxAPEC)]. The two 
plasma components represent typical high and low tem- 
peratures of a multi-temperature pla sma in an accretion 
column of an intermediate polar (e.g. Yuasa et al.[ 2010), 
and the Gaussian the fluorescence line from neutral Fe. 

In all cases, the fits successfully reproduced the ob- 
served spectra of these sources (Figure ^1), and pro- 
vided best-fit parameters as listed in Table [2] Model- 
predicted intensities in the 2 — 8 keV band were < 
2.4 x 10~ 12 erg cm" 2 s -1 (< 0.11 mCrab) for all the 
sources. When HXD/PIN observes these sources in its 



FOV, the angular transmission of the collimator (Taka- 



hashi et al. 2007) reduces these fluxes, and therefore their 
effective contributions become much smaller. 

Since typical NXB-subtracted HXD/PIN count rates 
in Regions 1 and 2 are ~ 0.04 and ~ 0.02 counts s _1 , 
respectively, a < 0.1-mCrab point source, with < 
0.003 counts s" 1 , contaminates the data at most a few 
to ten percent of the detected counts (before reduced by 
the collimator angular response). We consider that these 
sources are negligible in the HXD band compared to the 
NXB subtraction uncertainty (~ 0.003 counts s _1 at la). 

Other than these four sources, about ten point-source- 
like features are seen in the XIS images. Since their in- 
tensities, if treated as point sources, are all weak with 
fluxes on the order of 10~ 14 — 10~ 13 erg cm~ 2 s _1 in the 
2 — 8 keV band (equivalently 0.5 — 5 ^(Crab), their total 
flux amounts at most to a few percent of the detected XIS 
count rates. Therefore, we neglected these possible weak 
point sources when extracting the GRXE spectrum. 

Summarizing these evaluations, we conclude that the 
total unwanted contribution from the resolved point 
sources to the HXD/PIN data is considerably less than 
0.1 mCrab, or equivalently 3 x 10 -3 counts s^ 1 over the 
15 — 50 keV band. This means that the present analysis 
integrates X-ray signals emitted from unresolved sources 
fainter than 0.1 mCrab, or ~ 2 x 10~ 12 erg cm~ 2 s _1 in 
the 2 - 8 keV band, as the "GRXE" . This limiting flux 
corresponds to an intrinsic X-ray luminosity of a point 
source of 8 x 10 33 erg s _1 in the same energy band if a 
distance to the Galactic center of 8 kpc is assumed. This 
limit is severa l to ten times lower than those of the pre- 
vious studies ( |Revnivtsev et al.|2006| |Tiirler et aL|2010| ) 
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Galactic Longitude (deg) 



Fig. 1. — (a) A composite three-color X-ray image of the Galactic center and the Galactic bulge region created by adding data of 92 
mapping observations performed during 2005 — 2010. It was constructed from the XISO, 1, and 3 data. The 0.5 — 2, 2—5, and 5 — 8 keV data 
are represented by red, green, and blue, respectively. The image is corrected for the vignetting effect of the XRT and exposure differences, 
and the NXB was subtracted. Each r ectan gle represents a single observation. Labeled green crosses indicate resolved point sources of which 
the contributions were estimated in j ]2.4| (b) The same image as panel a but with the fill FOVs and ObsIDs of the utilized observations 
(Table [TJ overlaid; dashed and solid rectangles for the XIS and HXD/PIN, respectively. 
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in the hard X-ray energy band, thanks to the tightly- 
collimated FOV of the HXD and simultaneous imaging 
capability of the XIS. 

A further question might be raised about possible pres- 
ence of bright point sources which are inside some of the 
HXD/PIN FOV but outside the regions covered by the 
XIS. However, the INTEGRAL catalog already include 
dim point sources which were found in the ASCA Galac- 
tic center survey (Sakano et al. 2002), and we can con- 
clude that there is no other bright point source which 
could affect our analysis. 

2.5. Time Variability 

If any bright point source was inside the FOV of the 
selected observations, and if it varied on a time scale less 
than ^days, we should see count rate variations within 
individual observations. This issue applies particularly 
to the HXD, because it lacks imaging capability, and 
outer 75% (or 50%, after considering angular transmis- 
sion) of its FOV falls outside that of the XIS where no 
simultaneous imaging coverage is available. Therefore 
it is also important, before performing detailed spectral 
analysis, to confirm that the background-subtracted sig- 
nals of the HXD do not exhibit significant time varia- 
tions. This result also confirms that there was no bright 
fast transient point source which could contaminate the 
HXD data, i.e. located in the HXD FOV and outside the 
XIS FOV. 

We extracted light curves of HXD count rates as exem- 
plified in Figures p5l and examined them carefully. The 
periodic variation of the total count rates (ALL in the fig- 
ure) is mostly due to time variation of the NXB caused 
by changes of the cosmic-ray particle flux in orbit. After 
subtracting the NXB, residual count rates in these obser- 
vations (i.e. ALL-NXB) were ~ 0.02 - 0.07 counts s" 1 , 
and showed little time variation. Since the CXB amounts 
to ~ 0.01 counts s _1 in the HXD band, about 50 - 70% 
of these residual counts can be considered to originate 
from the GRXE. 

2.6. Extracting the GRXE Spectrum 

Following the inspection of data quality, we extracted 
spectra of the individual observations in Regions 1 and 
2. The XIS data of the entire imaging area of the CCD 
were used, except for two corners irradiated by calibra- 
tion radioactive isotopes and circularly-masked regions 
around the four point sources. The XIS NXB was esti- 
mated by applying the standard tool xisnxbgen to the 
same extracting region as the data. Although the XIS 
detects significant signals below 2 keV as well as above 
that energy, we did not used them because detailed spec- 
tral modeling in this energy band is hampered by com- 
plex interstellar absorption. As the HXD NXB, we used 
the simulated NXB file which is produced by the HXD 
team. Accuracies of these NXB subtractions are ~ 5% 



and 1% (la) in the XIS and the H XD, respectively (Tawa 
et al.||2008||Fukazawa et"al][2009| . Since the HXD NXFJ 
amounts to ~ 0.3 counts s in the 15 — 50 keV band, the 
uncertainty corresponds to ~ 0.003 counts s _1 . This is 
only 4 - 15% of the CXB-subtracted GRXE count rates 
in the same energy band (Table [TJ- 

Figure [4] shows examples of the spectra, derived from 
the same observations as were used in the previous light 



curve analysis. The raw spectra are shown as "ALL" , to- 
gether with "NXB" and "ALL-NXB" (NXB-subtracted 
celestial signals) spectra. As expected from the light 
curves, the NXB subtraction has left signals which are 
significant compared to the statistical and systematic un- 
certainties of the NXB subtraction. We regard this as the 
GRXE signal, again with 30 — 50% contribution from the 
CXB. A difference of spectral slopes below ~ 3 — 4 keV 
between the two regions is caused by the stronger inter- 
stellar absorption on the Galactic plane (Region 1) than 
that in the off-plane regions (Region 2). The intense Fe 
Ka emission lines (6 — 7 keV) are noticeable even in the 
individual XIS spectra. In addition, several emission- 
line like features are recognizable in 2 — 4 keV (e.g. at 
around 2.5 and 3.2 keV) which are thought to originate 
from lighter elements. 

Within each Region (1 or 2), the spectra of individual 
observations did not differ from one another within sta- 
tistical errors. Therefore, we combined their data for de- 
tailed broad-band spectral analyses, and created summed 
spectra of the two regions. The derived spectra are shown 
in Figure [5] The data summation has improved the 
data quality significantly, with little change in the spec- 
tral shapes. Based on the NXB-subtraction accuracies, 
we consider that the celestial signals are detected up to 
10 keV and 50 keV in the XIS and HXD/PIN, respec- 
tively. Figure [6] gives a close-up view of the Fe emission 
lines in the 6—7 keV band in Figure [5j 

In Region 2, we limited our spectral analysis to nar- 
rower energy ranges (2—9 keV for the XIS and 15—40 keV 
for HXD/PIN) due to lower counting statistics, which, in 
turn, is because of the shorter exposure and the lower sur- 
face brightness of the GRXE in off-plane regions. Com- 
pared to the individual spectra, the emission lines from 
Fe and the other lighter elements are more clearly seen 
in the XIS spectra. 

2.7. Surface Brightness of the GRXE 

To conduct spectral analyses, energy responses and ef- 
fective areas of the XIS and HXD /PIN should be calcu- 
lated for individual observations, by taking into account 
a global distribution of the GRXE surface brightness. 
This quantity is known to correlate with that of near 
infrared diffuse emission which reflects the stellar mass 
density ( Revnivtsev et al.|2006[ ) . Among many plausible 
func tions that can model the distribution, we also used 
after Revnivtsev et al.| ( |2006[ ), so-called G3 model (|Dwek 
et al. 1995) for the bulge/ bar structure. The model de 
scribes the stellar density p(r) as 



p(r)=p r 18 exp(-r 3 ), 



-1 1/2 



(1) 
(2) 



Where x, y and z are three-dimensional coordinates cen- 
tered on the Galactic center, and z axis corresponds 
to the pole of the Galactic rotation. See Dwek et al. 
( |1995| for detailed explanation of the assumed coordi- 
nate system, and a rotation of the bar structure. For 
three scale-length parameters, yq, ?/o, and zq, we used 
values presented in Dwek et al. (1995); xq = 4.01 kpc, 
yo = 1.67 kpc, and zq = 1.12 kpc. 
The surface brightness S at a certain sky position is 



(a) Soft point source at (l,b)=(-1.47, -0.33) 
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(b)AX J1742.6-3022 at (l,b)=( - 1.56,-0.20) 
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Energy (keV) 
(c) Suzaku J1740.5-3014 at (l,b)=( -1.70,0.25) 
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Energy (keV) 
(d)IGRJ17391-3021 at (l,b)=( -1.93,0.45) 




Energy (keV) 



Energy (keV) 



Fig. 2. — Background-subtracted and response-inclusive XIS spectra of the four contaminating sources recognized in Region 1; (a) a 
newly-found soft X-ray point source, (b) AX J1742.6-3022, (c) Suzaku J1740.5-3014, and (d) IGR J17391-3021. Crosses and solid curves 
in upper panels are data and the best-fit models, respectively. Lower panels present fitting residuals in terms of \- Inset labels represent 
model-predicted 2 — 8 keV energy fluxes from the sources, followed by the same value scaled by that of the Crab nebula. 
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assumed to follow a line-of-sight integral, 

/>oo 

S(l,b) = S p(x,y,z)ds, 
Jo 



(3) 



where (7, b) is a sky position in the Galactic coordinate. 
Sq and s are a normalization parameter and a line-of- 
sight distance to (a;, y, z) measured from the Sun, re- 
spectively. In this calculation, a distance between the 
Sun and the Galactic center was assumed to be 8.5 kpc. 

Based on this, we constructed a surface brightness 
map, and used it as an input in calculating the detector 
responses. To justify the procedure, we compared the 
actually measured HXD/PIN count rates of individual 
observations with the near-infrared surface brightness of 
Eqn. ([3| convolved with the HXD/PIN angular response. 
We chose HXD/PIN because the interstellar absorption 
is negligible unlike below 10 keV, and the CXB count rate 
can be accurately subtracted (whilst this is not possible 
in the XIS due to the absorption) . A nice correl ation was 
obtained as plot ted in Figure[7| Thus, we, after |Revnivt-| 
sev et al. (20061, reconfirmed that the usage of the model 



is appropriate m the response calculation. 

PARAMETERIZING THE HARD X-RAY SPECTRUM OF 
THE GRXE 
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TABLE 2 

The best-fit parameters for the contaminating point sources. 





Source 


Model a 


AT H 
cm 


r b 


keV 


X 2 /v 


-^2-8 kcV d 
mCrab 


a 
b 
c 

d 


New source 
AX J1742.6-3022 
Suzaku J1740.5-3014 

IGR J17391-3021 


PL 
PL 
2CIE 

PL 


2 4+ 3 ' 7 
13+j 

q 7+1.0 

^-1% 

4 o+ - 8 
4 - u -0.6 


3 1+ 15 
°- L -0.7 
9 7+O.3 

1 g+0-2 
1,M -0.1 


10.8(> 6.5) 
31(>4) 


1.20(25) 
0.91(28) 
0.81(86) 

1.32(25) 


0.01 
0.07 
0.07 

0.11 



a PL = an absorbed power-law model. 2CIE = two collisional ionization equilibrium plasma plus a 
gaussian model. 

b Photon index of the power-law model. 

c Plasma temperature. Numbers in brackets are lower limits on the 90% confidence level. 
d Model predicted flux in the 2 - 8 keV band (lCrab= 2.2 X 10 -8 erg cm' 2 s" 1 ). 
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Fig. 4. — Examples of broad-band X-ray spectra obtained in 
individual observations; Obs.IDs 503021010 and 502059010 repre- 
senting Regions 1 and 2, respectively. Black crosses are the celestial 
signal counts, which consist of the GRXE and the CXB, derived 
by subtracting the NXB (blue curves) from the raw counts (green 
crosses). Blue lab els specify atomic emission lines clearly seen in 
the XIS NXB (see |Tawa et al.|2008| . For clarity, spectra of XIS0 
are plotted, although the other two XIS cameras give fully consis- 
tent data. 



3.1. Fit with the IP Spectral Model 

Magnetic CVs, including IPs in particular, are thought 
to be a major contributor to the GRXE in the hard X- 
ray band because of their relatively low luminosi ties, high 
volume densities, and the hard spect ral shapes ( Revnivt 



sev et al. 2006 



20061. This has been di 



Sazonov et al 

rectly c onfirme d by Chandra in the Galactic center region 
( Muno et al.|2004 ). X-ray spectral shapes of IPs and Po 



lars are quite similar, both having the multi-temperature 
thermal nature; the only difference is that the plasma 
temperatures of Pola rs are lower due to e nhanced cy- 
clotron cooling (e.g. Cropper et al. 1998), the GRXE 
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Fig. 5. — Spectra of the GRXE (including the CXB) summed 
over all observations in Region 1 (top panel) and Region 2 (bottom 
panel), presented in the same way as Figure [4] 



spectra integrated above 15 keV should approximately 
be reproduced by the IP spectral model . 

To test this idea, we fitted the HXD /PIN spectra of Re- 
gions 1 and 2 with the IP spectral model which we numer- 
ically const ructed and verified in our previous study of 
nearby IPs ( Yuasa et al.|2010 ). This spectral model rep- 
resents multi-temperature plasma emission from an ac- 
cretion column on top of the magnetic poles of an IP. The 
WD mass and Fe abundance, which are its primary free 
parameters, determ ine the spectral shape. Compared to 
Yuasa et al. (2010), we updated the model by taking into 



account a plasma co oling function recently published by 
Schure et al.| ( |2009[ ) so as to ac curately mod el accreting 
gas with sub-solar abundances (Yuasa 2011). Although 
this affected little the spectral shape, it has increased the 
self-consistency of our model. 
In the fitting, the models for Regions 1 and 2 were 
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Fig. 6. — A close-up view of the same GRXE spectrum of Region 
1 as Figure[5] Three Ka lines from Fe in different ionization states 
are clearly resolved; Fe I = neutral, Fe XXV = He-like, and Fe 
XXVI = H-like. 
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Fig. 7.— Correlation of the HXD/PIN GRX E count rate 
(15 — 50 keV) and the near IR surface brightness (Dwek et al. 
|1995| convolved with the HXD/PIN angular transmission (34' X 34' 
at the full width at half maximum). Count rates of individual 
observations of Regions 1 and 2 are plotted in filled rectangles 
and open circles, respectively. Dashed line is the best-fit linear 
function (y = 4.6x — 0.005). The estimated CXB count rate of 
0.016 counts s — 1 was subtracted from the HXD/PIN count rates. 



constrained to share the same WD mass parameter, and 
allowed to have different normalizations. The Fe abun- 
dance parameter was fixed at unity in the fitting because 
it cannot be constrained without data of the emission 
lines. Even when we changed this parameter, for exam- 
ple, to 0.5 solar, the result was not affected at all. The 
CXB contribution is considered as a fixed model com- 
ponent using the esta blished models (e.g 



Boldt 1987 



Revnivtsev et al.|2003 ), rather than subtracting from the 
data, because the CXB signals are affected by interstellar 
absorption which is to be determined. The best-fit model 
successfully reproduced the overall spectral shape in the 
two regions, and gave xl = 1-37(30) with a null hypoth- 
esis probability of 9%. Therefore, we consider this fit as 
acceptable. Figure [8] shows the spectra, with the best-fit 
model superposed. The best-fit WD mass parameter is 
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Fig. 8.— The HXD/PIN spectra of the GRXE of Region 1 
(filled rectangles) and Region 2 (open circles), simultaneously fitted 
with the IP spectral model. Gray solid curves are the best-fit 
model spectra (CXB inclusive). Crosses in the lower panel show 
fit residuals. 



3.2. Fit with simpler models 

To characterize the hard X-ray GRXE spectra using 
simpler alternatives than the IP model, we also fitted the 
same HXD /PIN data with a CIE plasma emission model 
(the APEC model in XSPEC) and a power-law function. 
Although these fittings are rather empirical, the obtained 
representative spectral parameters are considered to be 
useful for comparison with previous studies in similar 
energy bands. 

A CIE plasma emission model yielded an acceptable fit 
with xt = 1-37(31), and the best-fit plasma temperature 
of kT — 15.77120 keV- If we simply take this value as 
a shock temperature, and convert it to a WD mass us- 
ing the relation between shock temperature and the WD 
mass in an IP ( Yuasa et al.| 2010| |Yuasa| 
mass of 0.48 M is derived. 



201 ip , a WD 
'his \VD mass is lower than 
that obtained by the IP model fit, because the "color 
temperature" of the IP model is generally lower than the 
actual shock temperature due to the multi-temperature 
nature of the post-shock region. 

A power-law fit to the HXD/PIN data gave a slightly 
worse but acceptable fit with xt = 1-43(31). The best-fit 
power-law index is T — 2.8 ± 0.2 which is slightly softer 
than the value d erived from the RXTE an d CGRO/OSSE 
data (2.3 ± 0.2; |Valinia & Marshall||1998| . If this power- 
law function is extrapolated down to the XIS energy 
band, the model-predicted XIS count rates exceed the 
actual data by more than a factor of 2 (i.e. too steep to 
reproduce the observed XIS spectra). Since this discrep- 
ancy will not be solved unless an artificial (and unrealis- 
tic) flattening of the power law in < 10 keV is assumed, 
we consider that the power-law fit is not an appropriate 
model to interpret the broad-band spectrum, although it 
roughly reproduces the HXD/PIN spectra. 

4. SPECTRAL ANALYSIS OF THE GRXE IN THE SOFT 
X-RAY BAND 
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TABLE 3 
Best-fit parameters for individual 
emission lines seen in the grxe 
spectrum. 



Line 


Energy 


cr a E.W. b 




(keV) 


(eV) (eV) 



S XV Ka 


2.45 


± 


0.01 


< 


35 


25+1 


S XVI Ka 


2.63 


± 


0.01 


< 


25 


35+1 


Ar XVII Ka 


3.13 


± 


0.01 


< 


11 


4ii 6 6 


Ar XVIII Ka 


3.33 


± 


0.02 


< 


10 


14+-; 


Ca XIX Ka 


3.86 


± 


0.07 


< 


89 


27+1 


Fe I Ka 


6.41 


± 


0.01 


< 


36 


76+ 9 
'°-10 


Fe XXV Ka 


6.68 


± 


0.01 


< 


17 




Fe XXVI Ka 


6.97 


± 


0.01 


< 


31 





1 Upper limits (90% confidence level) of in- 
trinsic line width in terms of 1 a of Gaussian. 

2 Equivalent width calculated against a con- 
tinuum containing the GRXE and the CXB 
signals. 

4.1. Identification of Atomic Emission Lines 

As can be seen in Figure [5j the XIS spectra ex- 
hibit many emission lines from astrophysically abundant 
heavy elements. To identify elements emitting these 
lines, as well as ionization states, we fitted the XIS 
spectrum with a phenomenological model consisting of 
a power-law continuum subject to an interstellar absorp- 
tion, and multiple Gaussians. As shown in Figure [9j the 
model successfully reproduced the spectrum when we in- 
troduced 8 Gaussians for prominent line features. The 
best fit center energies of these lines and their equivalent 
widths are listed in Table [3| With these center energies, 
identification is straightforward, and we found that He- 
like or H-like S, Ar, Ca, and Fe ions are emitting these 
lines as labeled in the figure. 

In the 2 — 3 keV energy range, the foreground diffuse 
X-ray emission is know n to contribute to some extent to 
the emissio n lines (e.g. Kuntz fe Snowden|2008 Yoshino 



et al. 2009 



Masui et al'|2009|) . Ho wever, the energy flux 



attributable to this emission is only a few percent of the 
detected flux of the GRXE, and therefore, is negligible 
at this stage. Instead, in a broad-band spectral analysis 
presented below, we take into account this foreground 
emission. 



spectra with a single-temperature CIE plasma model. 
The CXB contribution was included as a fixed model, like 
in the HXD/PIN analysis. As shown in Figure 10 (a), the 
best-fit model (kT ~ 4 keV) significantly underpredicts 
the S XV Ka (2.44 keV) and Ar XVII Ka (3.13 keV) line 
fluxes, and is unacceptable with \t — 2.91(245). These 
deficits cannot be resolved even when the abundance pa- 
rameter is increased, because the model would then over 
predict He-like lines of these elements. This reconfirms 
that the GRXE spectrum cannot be explained by the 
single-t emperature C I E pla sma model. 



After Kaneda et al. dl997|), we introduced an additional 
lower-temperature CIE component to better reproduce 
the low-energy lines. The fit improved considerably to 
X 2 = 1.46(243), and yielded the best-fit temperatures of 
fell = 0.79±g;|g keV and kT 2 = 6.23^;^ keV. Although 
this fit, with a null hypothesis probability less than 1%, 
is not yet acceptable in the strict sense, the remaining is- 
sue can be ascribed to calibration inaccuracy rather than 
to model inappropriateness, because the fit residuals are 
mainly observed in line- wing regions. If we include a 
model systematic error of 2% to simulate the insufficient 
response accuracy, the fit becomes formally acceptable 
with a null hypothesis probability of 1%. 

The derived two plasma temperatures are just repre- 
sentative ones, and may be approximating a mixture 
of more than two distinct plasma temperatures. The 
plasma may even have a continuous temperature distri- 
bution. However, the data statistically do not require 
any additional component. Therefore, instead of apply- 
ing more complex models to the XIS data, we proceed to 
broad-band spectral analysis combining the XIS and the 
HXD/PIN data. 

5. BROAD-BAND SPECTRAL DECOMPOSITION 



As shown in §3.1[ the hard X-ray spectrum of the 
GRXE is well reproduced with spectral models which 
have convex spectral curvatures, or more physically, the 
thermal nature (i.e. the IP model or a CIE plasma 
model). Besides, we reconfirmed, in the soft X-r ay band, 
the multi-temperature nature of the GRXE (Kaneda] 
et al. 1997|. Based on these, we try, in this section, to 



reconstruct the broad-band GRXE spectrum (Figure [5| 
using as small a number of physically plausible spectral 
components as possible. 



4.2. Multi-temperature Nature of the GRXE 

It is almost obvious from the above identification that 
at least two distinct plasma components compose the 
GRXE. This is because, as long as collisional ionization 
equilibrium is assumed, elements like S and Ar would be 
almost fully ionized and would not emit emission lines 
if the plasma is hot enough (10 7 5-8 K) to highly ion- 
ize Fe up to Fe XXVI Ka (6.9 keV; H-lik e); see ioniza- 
tion b alanc e calculated by , for e xample, |Bryans et al.| 
( |2006[ ) and |Bryans et al.| {2009). This reconfirms the 
multi-tempera t ure na ture of the GRXE first revealed by 
Kaneda et al. (1997) based on similar arguments using 
AbUA/GiH data. In addition, the existence of the Fe 
I (neutral) Ka implies reprocessing of X-rays by a cold 
matter. 

To quantify the multi-temperature nature of the 
GRXE found above, we first performed a fit to the same 



5.1. Model Construction 

In analyzing the 2 — 50 keV GRXE spectrum, let us 
choose, as our starting point, the two CIE plasma com- 
ponents which gave a successful fit to the XIS spec- 
tra. In addition to them, as noted above, we con- 
sider a contribution from the foreground diffuse soft X- 
ray e mission. Based on r ecent rep orts ( |Yoshino et al 



2009 Masui et al. 2009 Kimura 2010), we assumed 



the surface brightness of this foreground emission to be 
2.5 x 10 -9 erg cm -2 s _1 sr _1 , and included it in the 
spectral model as a CIE plasma component. The tem- 
perature was fixed at kT = 0.7 keV, because the fore- 
ground emission contributes no more than 1%. Here- 
after, we designate a sum of the two CIE plasma emission 
plus the foreground emission as "Model 1" for simplic- 
ity. Since the two plasma components in the GRXE can 
come through different absorption column densities, they 
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Fig. 9. — Emission lines seen in the GRXE spectrum of Region 1 taken with XISO (black crosses). Plotted model (gray solid curve) is a 
phenomenological one consisting of a power-law continuum (black dotted curve) and eight Gaussians (black solid curves) for emission lines. 
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Fig. 10.— (a) The GRXE spectrum of XISO (black crosses) fitted 
with a single-temperature CIE plasma emission model and addi- 
tional components, i.e. the CXB (gray dashed curve) and a Gaus- 
sian for the Fe I Ka line (gray dotted curve). Red and magenta 
curves are the total model and the CIE plasma component, respec- 
tively. Vertical dotted lines indicate those emission lines which 
are not reproduced by the model. Fitting residuals are shown in 
the lower panel, (b) An improved fit incorporating two single- 
temperature CIE plasma emission components (magenta and blue 
curves). The CXB in this panel has a spectral shape which is 
different from that of panel a below 3 keV because the best-fit in- 
terstellar absorption column densities differ (JVh 
in panel a, while _/V H = 4.8 X 10 22 cm -2 in b). 



x (Foreground + Abs.l x CXB 
+Abs.2 x CIE1 + Abs.3 x CIE2). 

The 6.4-keV Fe Ka emission line could provide a pow- 
erful tool for understanding the origin of the GRXE be- 
cause, near a WD, its line shape is probably distorted 
by Compton scattering in the atmosphere or the gravi- 
tational redshift. However, these studies would require 
a higher energy resolution. In the present analysis, we 
therefore approximate it by a Gaussian with its normal- 
ization left to freely vary. 

Another model, Model 2, was constructed by replac- 
ing one of the two CIE plasma com ponents of Model 1 
with our IP spectral model (IP_PSR; |Yuasa et ai]|2010[ ). 
This can be supported by the successful appl icati on of 
this model to the HXD/PIN GRXE spectrum fljO), and 
also by a suggestion that magnetic CVs are a major con 
tributor to the GRXE above 10 keV (e.g. Revnivtsev 



et al. 20061. The Fe abundance of the IP component 
was tied with that of the CIE Plasma component. Other 
abundance parameters were again allowed to individually 
vary (like Model 1). Since intrinsic par tially-covering ab 



sorpt i on is a common f eature of IPs ( Suleimanov et al. 
2005| Yuasa et al.||2010 ) , the IP model was subjected to 
an additional (partially-covering) dense absorption, with 
a covering fraction /, of which absorption column den- 
sity was allowed to freely vary. Summarizing above, this 
Model 2 is constructed as 

Model 2 = Abs.l x (Foreground + Abs.l x CXB 
+ Fe I Ka + Abs.2 x CIE 
+ [/ ■ Abs.3 + (1 - /) ■ Abs.4] x IP_PSR). 

We further subdivided this model into Model 2a and 
2b, by treating the AfwD parameter of the IP component 
differently. In Model 2a, Mwtj was fixed at 0.66 M de- 



rived from the HXD/PIN spectral fitting (£ 3.1 1 . In con- 
trast, the WD mass parameter of Model 2b was allowed 
to freely vary. 



were subjected to independent absorption factors. The 
abundances of major elements (i.e. Si, S, Ar, Ca, and Fe) 
were allowed to vary individually, but were constrained 
to be the same between the two CIE components. Thus, 
Model 1 can be expressed as, 

Model 1 = Abs.l (4) 



5.2. Combined Spectral Fitting 

We first fit the models to the stacked GRXE spectrum 
of Region 1, because it has higher counting statistics and 
covers a wide r energy range over 2 — 50 keV. As drawn 
in Figure [IT] all the three models well reproduced the 
broad-bandspectrum, and gave the best-fit parameters 
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as expla ined in §4.2 
errors (jSjj MoHe 



as listed in Tables [|] and [5j Although the values of x 2 
are not yet small enough to make the fits formally ac- 
ceptable, we consider again that this is caused by the 
inaccuracy of t he X IS response in the 2 — 4 keV band 
If wc include additional systematic 
2 became acceptable with null hy- 
pothesis probabilities larger than 1%, while Model 1 was 
still unacceptable (with a probability of 0.2%). 

The CIE1 component of Model 1 and CIE of Model 
2 yielded a plasma temperature of kT — 1.4 — 1.7 keV, 
and account for the soft continua plus the Fe XXV Ka 
line with (almost) no contribution to the Fe XXVI Ka 
photons. Most of the hard X-ray flux detected in the 
HXD/PIN band is explained by the higher temperature 
components, CIE2 with a temperature of 15.1^Q'y keV 
in Model 1, or the IP component with WD masses of 



j+0.05 
-0.04 



Mq in Model 2a and 2b, respec- 



0.66 M Q and 0.48Z 
tively. 

Similarly, the stacked Region 2 spectrum was fitted 
with the same model. Due to poorer data quality, we 
fixed some parameters which mainly determines spectral 
shapes of the individual components; the CIE plasma 
temperatures and Mwd as summarized in Tabl es [4] and 
[5j The best-fit models are plotted in Figure T2| 

Figure [13] shows the best-fit Model 2b in Region 1 by 
removing the foreground emission, the CXB, and the in- 
terstellar absorption. A hump-like structure seen in the 
10 — 30 keV band is due to the partially-covering dense 
absorption applied to the IP model. 

6. DISCUSSION 

6.1. Interpreting the Obtained Spectral Parameters 

As presented in the previous section, the GRXE spec- 
trum up to 50 keV can be well described only with ther- 
mal plasma models without any empirical power-law con- 
tinuum. We consider that this is a strong support for our 
working hypothesis, the "Point Source" scenario of the 
GRXE origin because the derived plasma temperatures 
and WD masses are well consistent with values expected 
from this scenario. Contrarily, in the context of the "Dif- 
fuse" scenario, problems of confinement and energy in- 
jection must be solved for sustain hot interstellar plasma 
with a temperature higher than 10 keV in the Galactic 
plane. 

The l owe r-te mpe rature CIE component (magenta in 
Figures 11 and 12) resulted in ranges of kT — 1.4 — 
1.7 keV and 1.2-"L3 keV in Region 1 and 2, respectively. 
These temperatures are well consistent with tho se seen 
in X-ray spectra of coronal X -ray sources (e.g. 



Favata 



et al.|1997 Covino et al.|2000| ) that are binary stars har 
boring two late-type (or "normal" ) stars like the Sun, 
and exhibits X-ray emission via magnetic activities in 
their coronae. The small temperature difference between 
the two regions (~ 0.3 — 0.4 keV) can be naturally un- 
derstood if there is unresolved contribution from young 
supernova remnants (SNRs) which tend to exhibit ther- 
mal X-ray spectra of pl asma (electron) temperatures of 



kT ~ 2 - 4 keV (e.g. Kinugasa fe Tsunemi 1999 for 



Kepler's SNR; Tamagawa et al.||2009 for Tycho's SNR). 



99] 
SNI 



Recent discoveries or possible discove ries of new SNRs in 
the Galacti c center region ( Sawada et al. 2009 Nobukawa 
et al. 2008 1 , especially in on-plane regions such as Region 
1, support this speculation on the contribution from un- 
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Fig. 11.— The GRXE spectrum of Region 1, fitted with Model 
1 (top panel), Model 2a (middle panel), and Model 2b (bottom 
panel). Black crosses are the observed spectrum (NXB subtracted), 
and red curves are the sum of all model components. In the Model 1 
fit, the two CIE plasma components are shown in magenta (lower 
temperature) and green (higher temperature). Magenta and or- 
ange curves in the lower two panels are the CIE plasma and the 
IP model components, respectively. Solid, dashed, and dotted gray 
curves are components representing the foreground diffuse soft X- 
ray emission, the CXB, and Ka emission line from neutral Fe, 
respectively. 



resolved SNRs. 

The temperature of the CIE 2 component of Model 1, 
15.llg' 7 keV, is consistent with a representative plasma 
temperatu re (10 — 20 keV) of magnetic CVs (e.g. |Yuasa| 
et al. 2010| ) . The Model 2a fit successfully reproduced the 
broad- band spectrum although a WD mass parameter 
was fixed at the val ue de rived from the HXD/PIN anal- 
ysis (0.66t u ;^ M Q ; ^3T|. When it was unfixed in Model 
2b, a similar but slightly lower WD mass of 0. 48lg; U 4 M 
was obtained. This reduction in WD mass is also ex- 



plained by unresolved contributor from other types of 
CVs, especially dwarf novae in quiescent, whose spectra 
can also be reproduced with a CIE plasma model with 
lower representative temperatures (kT ~a few— 10 keV) 
than those of magnetic CVs. 

Is the reduction of the WD mass parameter by ~ 
0.10 — 0.15 Mq realistically possible with the contribu- 
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TABLE 4 

Result of the wide-band spectral fitting with Model 1 





n H a 


fcT CIE1 b 


n H c 


fc7ciE2 b 


7 cl 




2 c 
*k, 3% 


^CIEl 


^CIE2 




10 22 


cm 


keV 


10 22 cm~ 2 


keV 


Zq 






(2 - 10) 


(2 - 50) 


Region 1 


4.0 


+0.2 
-0.2 


1 6fi+ - 04 


I5.9ti; 5 3 


K 1+0-4 
io ' i -0.7 


n 07+0.06 

u - y '-0.06 


1.45(754) 


1.15 


1.9 


7.8 


Region 2 


1.3 


+0.2 
-0.2 


i 01 +0.03 
1 - J± -0.03 






i ofi+o- 07 

4-UD_ og 


1.85(353) 


1.25 


7.9 


11.8 



Blank parameters are fixed at the Region 1 values (see text). 
a Hydrogen column density of interstellar absorption. 
b Plasma temperatures of the CIE 1 and the CIE 2 components. 
c Additional absorption column density applied to the CIE 2 component. 
d Fe abundance of the CIE components. 

c Improved fitting statistics achieved when the XIS response uncertainty was virtually took into account. 

f Model-predicted fluxes of the CIE 1 and 2 components in 10 -9 erg cm -2 s _1 integrated over specified energy ranges and 

a \l\ < 10° and |6| < 10° region. 



TABLE 5 

Result of the wide-band spectral fitting with Models 2a and 2b* 





n H a 
10 22 cm -2 


kT 
keV 


Z Fe h 

Zq 


10 22 cm" 2 


M WD d 




2 e 
K v, 3% 


^CIE f 
(2 - 10) 


(2 - 50) 


Model 2a 
Region 1 
Region 2 
Model 2b 
Region 1 
Region 2 


3 6+0-2 
i O+0.2 

9+0.4 

-0.6 

1 l+°- 2 


-, c-9+0.04 
1 - oz -0.04 
i 91 +0.04 
1 - zl -0.03 

1 44+0-06 
1 - 44 -0.07 
-, i 7 +0.04 
1 - 1 ' -0.04 


o.8«t8: o$ 

0.82t°;°° 

n 70+O.O6 

67+0' 04 
u - D '-0.04 


13.0l?J 
ll.8±f;| 


0-66( flxed ) 

n 4s+0- 05 
U.4S_ 04 


1.30(753) 
1.62(354) 

1.28(752) 
1.65(354) 


1.04 
1.11 

1.02 
1.13 


3.9 
7.2 

3.2 

6.5 


8.6 
12.7 

8.6 
12.5 



Blank parameters are fixed at the Region 1 values (see text). 
a Hydrogen column density of interstellar absorption 
b Fe abundance. 

c Additional absorption column density applied to the IP component 
d WD mass of the IP model. 

° Improved fitting statistics achieved when the XIS response uncertainty was virtually took into account. 

f Model-predicted fluxes of the CIE and IP components in 10 — 9 erg cm~ 2 s" 1 integrated over specified energy ranges and a 
\l\ < 10° and |6| < 10° region. 



tion of dwarf novae? To test this, we simulated a com- 
posite broad-band GRXE spectrum (excluding the low- 
temperature CIE component) by adding spectral models 
of dwarf novae and IPs. A representative plasma tem- 
perature of kT = 5 keV and an average WD mass of 
0.6 Mq were assumed for the dwarf nova and the mag- 
netic CV components, respectively. Relative fluxes of 
individual components were adjusted according to X-ray 
emissiyities of these source types measured by |Sazonov| 
et al. (2006). For better simulating the observed spec- 
trum, we also applied photo absorption models with hy- 
drogen column densities of Ah = 3 x 10 22 cm~ 2 (Galac- 
tic interstellar value for dwarf novae) and 10 x 10 22 cm -2 
(dense intrinsic absorption for magnetic CVs). Simulated 
count rates, or absolute intensities, were a dju sted so as to 
match those of the observed data. Figure [14] (a) presents 
a thus produced spectrum. This spectrum mimics the 
GRXE spectrum (subtracted the low-temperature CIE 
component). 

We performed a similar fit to this simulated spectrum 
using only the IP model, and obtained an acceptable fit 
which yields xL= 1.13(287). The best-fit model is plot- 
ted in Figure 14 (b), and it gave the best-fitting WD 
mass of 0.51 ± -0701 Mq and absorption column density 
of A H = 7.6 ± 0.3 x 10 22 cur 2 . Thus, even when the 
dwarf nova component is missing in the fit model, the 
composite simulated spectrum can be well reproduced 



with the IP model with the WD mass parameter which 
is slightly reduced from the assumed value of 0.6 M Q 
(AMwd = 0.09 M Q ). This supports the above discus- 
sion on the reduced WD mass obtained in the broad-band 
spectral analysis. 

6.2. The Number Density of the Unresolved Hard X-ray 
Point Sources 

The GRXE surface brightness can be used to con- 
strain population of unr esolved point sources, mostly IPs 
( Revnivtsev et al.|20 06). Using the GRXE flux observed 
by HXD/P1JN , we derive their population, or so-called X- 
ray luminosity function. In the present calculation, we 
assume that all the point sources are located at a dis- 
tance of 8 kpc from the Sun (i.e. at the Galactic center 
region) as done in many previous studies (e.g. Muno 
|et al.||2009[ ). 

Generally, a luminosity function of X-ray stars A can 
be expressed as a power-law function, and therefore we 
write A using the intrinsic luminosity L 



N{> L) = N 



(5) 



Here, A(> L) means the surface number density of (un- 
resolved) point sources which have an intrinsic luminosity 
L, and has a dimension of [(number) (solid angle)" 2 ]. Ao 
and Lq are constant factors for scaling, and the latter can 
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U 1 



Model 1 : CIE_Plasmal + CIE_Plasma2 



Region 2 





Model 2a : CIE_Plasma + IP_PSR (M wu =0.66 M sun fixed) 



Region 2 





Model 2b : CIE_Plasma + IP_PSR (M WD free) 



Region 2 






(a) Simulated spectrum of dwarf novae and magnetic CVs 



10 

Energy (keV) 



20 



50 



Fig. 12. — The same as Figure fLll but for Re gion 2. The tem- 
arature U1F 



pcrature of the lower-temperat 
freely vary. 



IE component was allowed to 



+ Krivonos et al. 2007 (RXTE/PCA) 

+ Krivonos et al. 2007 (INTEGRAL/IBIS) 

+ Ttirler et al. 2007 (INTEGRAL/IBIS) 




10 20 

Energy (keV) 



200 



Fi g. 13. — (a) The same best-fit GRXE model spectrum as Fig- 
ure [Tl] but represented by removing the detector responses. Data 
points take n from previous stud ies are overlaid; green a nd blue 
crosses from lKrivonos et al.|f2007) and magenta crosses fromlTiirlerl 

Etfqp5i(f 11 — ' 



be arbitrarily fixed. The parameter a denotes luminosity 
dependence of N. 

Now, we can calculate an X-ray surface brightness S 
predicted from the source population N by integrating a 
product of the number of sources (i.e. derivative of the 




magnetic CVs 



dwarf r 



2 5 
(b) Fit with the IP model 



10 

Energy (keV) 




n H -7.6(7.3-7.9) 10 22 cm" 2 
M^-o.siOi.so-o.sz) M si 

X', =1.13(287) 





2 5 10 20 50 

Energy (keV) 

Fig. 14. — (a) A simulated GRXE spectrum (black crosses) 
composed of spectral models of the dwarf nova (blue curve) and 
the magnetic CV (orange curve; the IP model), (b) The same 
spectrum as panel a, but fitted with the IP model (orange curve). 



cumulative density N) and the luminosity L as 



S 



AL 



L dL 



(6) 



The integration is performed over a luminosity range 
^min — L max where dN(> L)/dL is non-zero. 

For quantifying the luminosity function from the 
present GRXE flux, we set £ m in an d i m ax as follows. 
The present GRXE study does not contain bright known 
X-ray sources which have energy fluxes higher than 
~ 1CP 12 erg cnr 2 s _1 in the 2 — 8 keV band, or 
equivalently ~ 0.1 mCrab. The flux limit corresponds 
to 1.7 x 10 -12 erg cm -2 s _1 when expressed in the 
15 — 50 keV band (i.e. the HXD/PIN energy cover- 
age) . If we consider that a point source located at 8 kpc 
from the Sun, the flux means an intrinsic luminosity of 
1.3 x 10 34 erg s _1 in the same energy range. This should 
be regarded as L max . As for L m - m , recent deep Chandra 
observations revealed that there exist X-ray sources with 
lumin osities as low as 10 30 er g s _1 in the 0.5 — 7 keV 
band (Revnivt sev et al.||2009[ ). Another important fact 
is that Muno et al.| ( |2U09| ) securely measured the shape 
of the luminosity function down to 3 x 10 31 erg s _1 in the 
2—10 keV band. This value provides stringent "upper 
limit" for L min ; i.e. its actual value should be much lower 
probably one order of magnitude because there is no sign 
of break, or turn off, in their luminosity function. In the 
present studywe tentatively considered two cases with 
i'min being set these measured values when integrating 
of Eqn. (61). If true L min is lower, the normalization of 
the calculated luminosity function will decrease. 

Before the integration, these luminosity values should 
be converted to those in the 15 — 50 keV band which 
we concentrate on in the present calculation. Since 
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most of the unresolved point sources in this energy 
band are thought to be accreting WDs (especially mag- 
netic ones), th ey have hard spect ra (i.e. high plasma 



temperatures; Muno et al. 2004). We assumed, for 



the unresolved sources, a typical X-ray spectrum con 
sisting of a single-temperature CIE plasma model with 
kT = 20 keV suffered from the interstellar absorption 
of iV H = 6 x 10 22 cm~ 2 (typical value for the Galac- 
tic center region). Based on this spectral shape, we 



obtain converted L m i n values of 2 x 10 30 
3.8 x 10 31 erg s" 1 in the 15 - 50 keV. 



erg s 



and 




10" 



This work 

Chandra Galactic center 
(Muno et al. 2009) 




l(f 1(T 10" 10" 10 J 

Luminosity (erg s^in the 12-50 keV) 

Fig. 15. — X-ray luminosity functions of unresolved point 
sources, which compose the GRXE in the higher energy band 
above 10 keV, calculated from the HXD/PIN GRXE flux assum- 
ing L m i n = 2 X 10 30 (green region) and 3.8 X 10 31 erg s — 1 (orange 
region) in the 15 — 50 keV band. Two solid lines which enclose 
each region have differently-assumed power- law indices a = 1.0 
and 1.5. Blue dashed line represents a luminosity function actually 
measured by Chandra in similar sky regi on (2° X 0.8° arou nd the 
Galactic center) in the 2-10 keV band | |Muno et al.|2009} . 



The present HXD/PIN measurement gives the 
GRXE surface brightness, for example, S = 2.34 x 
10 35 erg s" 1 (PIN FOV)- 1 in Obs.ID 504001010 cen- 
tered at (I, b) = (-1.47°, -0.26°) in the 15-50 keV band. 
Since S, or N(> L) inside it, includes two unknown pa- 
rameters a and Nq, an additional constraint on one of the 
two parameters s hould be placed to determine N(> L) . 
Previous studies ( |Ebisawa et al.|2005| |Muno et al.|2009 1 
have revealed that the luminosity function of faint Galac- 
tic X-ray point sources has a = 1.0 — 1.5. Based on this, 
we calculate N(> L) using assumed indices a = 1.0 and 
1.5. 

Figure [15] shows thus calculated luminosity functions 
of the unresolved sources. Two filled regions repre- 
sent two representative cases with lower luminosity lim- 
its of Lmin = 2 x 10 30 erg s _1 (green region) and 
3.8 x 10 31 erg s _1 (orange). The luminosity functions 



are well consistent with that actually measured in sim- 
ilar sky regions by Chandra (Muno et al. 2009). From 



this comparison, we consider that the number density 
of required faint point sources is not unrealistically high 
unlike what prop osed by some previous studies (e.g. Ebi 
sawa et al.||2005| ). 



6.3. The Origin of the GRXE 

The broad-band spectral decomposition has been an- 
ticipated for long time for understanding the emission 
mechanism(s) and the energy supplier(s) of the GRXE. 
This has been done in the present study, decompos- 
ing the GRXE i nto two representative constituents. As 
was detailed in £5.2 the low-temperature and the high- 
temperature CIE plasma emissions have plasma temper- 
atures of ~ 1 keV and > 10 keV according to our model- 
ing. The latter component can be successfully replaced 
by the IP spectral model as examined in the Model 2 
fits. In addition, the luminosity function of unresolved 
hard X-ray sources does not largely contradict with one 
directly measured with Chandra. Based on these, we 
consider that the present results in < 10 keV provide 
yet another very strong support for the "Point Source" 
scenario of the GRXE along s ide of the imaging decom- 



position of Revnivtsev et al. ( 2009 ) 



On the other hand, the results are against for several 
previous studies that suggested the "Diffuse" scenario as 
the GRXE origin. The putative non-equilibrium ioni za- 
tion (NEI) plasma proposed by Kaneda et al. ( 1997 ) is 
rejected since we detected intense Fe XXVI Ka line in 
the present sky region which is not expected from their 
NEI model. This is a confirmatio n of a similar denial 
argument by Ebis awa et al. |2008[) in t he di fferent sky 
region. Unlike " Yamasaki et al. 1 1997 ) and Valinia & 



Marshall (1998), the observed broad- band spectrum did 



not require a putative hard-tail component which was 
suggested to smoothly connect to the gamma-ray back- 
ground emission up to hundreds of MeV. 

6.4. The Mean WD Mass in the Galaxy 

The WD mass of CVs is an important parameter when 
interpreting the hard X-ray GRXE spe ctral shape. It 
was first estimated to be ~ 0.5 M m in Krivonos et al 



( 2007[ ) by ro ughly fitting their GRXE spectra with the 
IP model by Suleimanov et al. ( |2005[ ). Our hard X-ray 
analysis confirmed this, an d gave a slightly heavie r WD 
mass of 0.66 M n . As 



Krivonos et al. 



(20071 men- 



tioned, the value could be interpreted as the mean WD 
mass of cumulated magnetic CVs in the Galaxy. 

If we accept this idea, it is intriguing to compare the 
derived WD mass with that reported for isolated WDs 
based on op tical spectroscopy. F or example, using large 



SDSS data, Kepler et al 
mass of 0.593 ± 0.016 M e 



(2007) reported a mean WD 
for 1733 WDs which have he- 



lium and hydrogen outer layers. Our value 0.66 



+0.09 
-0.07 



includes the average within errors. The slight difference 
of the center values (~ 0.06 M ) could be a result of 
long-lasting mass accretion. However, CVs are thought 
to evolve through the common-envelope phase of main- 
sequence stars, in which two stars in a binary share their 
outer layers, possibly resulting in WD masses different 
from those of isolated stars. Therefore, at this moment, 
it is difficult to investigate the difference of the average 
WD masses. More accurate WD mass determinations 
in CVs, and sophisticated evolution model of them are 
necessary to accurately perform above comparison. 

6.5. Connection to the Galactic Center X-ray Emission 

A low-temperature CIE component (kT ~ 1 keV) sim- 
ilar to one seen in the GRXE was also reported recently 
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by Nobukawa et al. (2010) in the Galactic center X- 
ray emission observed in [t| < 0.2°. Probably due to 
lack of the hard X-ray spectral coverage in their anal- 
ysis, the authors reported a plasma temperature of the 
hotter CIE component of 7.0 ±0.1 keV, which is consid- 



erably cooler than our result kT = 15. keV. Note 
also that their spectral analysis introduced, in addition 
to two CIE plasma components, an intense power-law 
component which was not detected in our broad-band 
GRXE data. 

In the Galactic center region, the super-massive black 
hole (~ 3x 10 6 Mq) is suspected to have a close con- 
nection with the unresolved (or literally diffuse) emi ssion 
which is observed in its proximity (e.g. predicted b y Sun- 
yaev et al. 1993| and observed e.g. by jlnui et al. 2009 U 
This emission could have another origin which differ from 
that of the GRXE. Besides, stellar density is higher by 
more than an order of magnitude in the very central re- 
gion of the Galaxy (~ 100 pc from the Galactic center) 
compared to those of the region we studied. Therefore, 
X-rays from faint point sources (mostly coronal X-ray 
sources and IPs) might occupy a considerable fraction 
of the detected signals, although Nobukawa et al. ( 2010[ ) 
simply neglected them. Since about a half of such faint 
point sources radiate its energy in the hard X-ray band, 
we stress that a broad-band spectral analysis including 
the hard X-ray energy range is essentially powerful to ex- 
amine possible existence of truly diffuse emission in the 
Galactic center. 



We accumulated the data of Suzaku GRXE observa- 
tions in the Galac tic b ulge region achieving ~ 1 Ms 
exposure in total ( |2.3[ ), and produced the broad-band 
GRXE spectra with ' high counting statistics covering the 
2-50 keV band (pjj). We showed that the hard X- 
ray GRXE spectrum taken with H XD/ PIN is well repro- 
duced by the IP spectral model (£ 3.1 ) with a WD mass 



parameter of 0.66^Q Qy M Q . From emission line analyses 
based on the CIE modeling w hich is more approp riate 



than previous line studies (e.g. Kaneda et al 



1997), we 



7. SUMMARY 



r econ firmed the multi-temperature nature of the GRXE 

Based on the se re sults, we constructed physical models 
of the GRXE Q5.2) which consists of multi-temperature 
CIE plasma emissions. The models nicely decomposed 
the broad-band GRXE spectral shape for the first time. 
Especially, the model that includes the IP component, 
namely Model 2, gave better fits to the data together 
with the low-temperature CIE plasma component. The 
derived plasma temperature and the WD mass are quite 
consistent with those of coron al X- ray sources and typical 
of magnetic accreting WDs ({ 6.1 ). 

We also calculated the X-ray luminosity function and 
the number density o f the unresolved hard X-ray point 
sources, mostly IPs (^6.2). The luminosity function is 
consistent with that obtained from a deep imaging obser- 
vation near our field, indicating that no unknown X-ray 
source type is additionally required to explain the GRXE 
flux. 

Combining the spectral decomposition result and the 
calculated source density, we concluded that the present 
result supports the "Point Source" scenario of the GRXE 
origin. 



APPENDIX 

SUZAKU OBSERVATIONS IN THE GALACTIC CENTER REGION 

For later reference, we summarize information of 92 Suzaku observations in the Galactic center region in Table [6| 
XIS data of these observations are used to produce the mosaic image shown in Figure [T] The total effective exposures 
of the XIS and HXD/PIN are 4.4 and 3.9 Ms. The list contains observations with various types of scientific aims; the 
Galactic Ridge X-ray Emission, supernova remnants, transient X-ray binaries, molecular clouds reflecting X-rays, and 
unidentified sources in TeV wavelength (H.E.S.S unID sources). 
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TABLE 6 

Suzaku observations of the Galactic center. 



Obs. ID a Coordinate 15 Start time Exposure 







I 


I) 




UT 


XI s 


PIN 


1 


100027010 


0.06 


-0.08 


2005-09- 


-23 


07:18:25 


44.8 


37.9 


2 


100027020 


-0.24 


-0.05 


2005-09- 


-2d 


14:17:17 


42.8 


36.1 


3 


100027030 


-0.44 


-0.39 


2005-09- 


■24 


"i 1 n'7 no 

ll:U7:0o 


2.1 


1.9 


4 


100027040 


-0.44 


-0.07 


2005-09- 


-24 


12:41:33 


1.9 


1.8 


5 


100027050 


0.33 


0.01 


2005-09- 


-25 


17:29:12 


2.0 


1.8 


6 


100037010 


-0.24 


-0.05 


2005-09- 


-29 


04:35:41 


43.7 


39.4 


7 


100037020 


-0.44 


-0.39 


2005-09- 


-30 


r\ a on A A 

U4:3U:44 


3.3 


3.1 


8 


100037030 


-0.45 


-0.07 


2005-09- 


-30 


r\i~- r\r- on 

06:06:32 


3.0 


2.8 


9 


100037040 


0.06 


-0.08 


2005-09- 


■30 


n""7 i o n 1 

07:43:01 


43.0 


39.5 


10 


100037050 


0.33 


0.01 


2005-10- 


■01 


t \ t • r>n 1 1 

06:22:41 


2.4 


2.2 


11 


100037060 


0.64 


-0.10 


2005-10- 


•10 


■in no n 1 

12:28:01 


76.6 


70.8 


12 


100037070 


1.00 


-0.10 


2005-10- 


■12 


07:10:24 


9.2 


9.5 


13 


100048010 


0.06 


-0.08 


2006-09- 


-08 


nn no n a 

02:23:24 


63.0 


60.3 


14 


102013010 


0.06 


-0.08 


2007-09- 


-03 


19:01:10 


51.4 


44.5 


15 


402066010 


-1.93 


0.45 


2008-02- 


■22 


11:52:49 


36.5 


31.3 


16 


403001010 


1.36 


1.05 


2009-02- 


-22 


"i n n a "in 

19:04:19 


71.5 


59.7 


17 


403009010 


0.17 


0.03 


2009-03- 


-21 


nn no no 

02:03:28 


110.8 


91.7 


18 


500005010 


0.43 


-0.11 


2006-03- 


-27 


no nn nn 

z3:UU:2z 


88.4 


64.6 


19 


500018010 


-0.57 


-0.09 


2006-02- 
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TABLE 6 
Continued. 
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Coordinate 13 
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a Observation ID. 

b Aim point in Galactic coordinate (degree). 
c Net exposure in units of 10 3 s. 
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